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Preface 

This  report  is  the  sixth  concerned  with  research  accomplished  in 
connection  with  Navy  Contract  Nonr-433(00) , between  Dunham  Laboratory, 
Yale  University,  and  the  Office  of  Naval  Research,  Department  of  the 
Navy.,  In  this  report  is  given  a discussion  of  the  design  of  a class 
of  networks  that  can  be  used  for  producing  a known  constant  time  delay 
with  an  electronic  analog  computer,  A differential-difference  equa- 
tion of  the  sort  considered  in  Report  No,  5 requires  such  a delay 
network  if  it  is  to  be  studied  simply  on  a computer.  It  was  in 
connection  with  the  investigation  of  this  equation  that  a study  of 
delay  networks  was  made. 

The  research  was  carried  on  and  the  report  written  by  the 
undersigned. 

W,  «T»  Cunningham 
New  Haven,  July  1954 
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Abstract 


In  making  use  of  an  analog  computer  to  study  differential- 
difference  equations  in  which  a constant  time  delay  occurs,  it  is 
necessary  to  have  soma  device  to  introduce  a known  constant  delay  into 
the  system.  Usually  it  is  convenient  to  do  this  with  a network 
synthesized  from  the  elements  of  the  computer.  A perfect  delay  netvrork 
is  an  all-pass  network  in  which  phase  shift  varies  linearly  vrith 
frequency.  One  class  of  networks  that  approaches  thi3  ideal  is  studied 
by  considering  the  roots  of  their  transfer  functions,  Networks  having 
one,  two,  three,  and  four  pairs  of  roots  are  investigated  in  detail, 
and  numerical  data  about  them  nm  nw rented  in  tabular  for m*  These 
tables  allow  selection  of  a particular  network  appropriate  to  the 
problem  being  analyzed. 


1 


I.  Need  for  Delay  Networks 

Certain  physical  systems  have  within  them  some  process  that  intro- 
uubos  a time  uelay  of  constant  magnitude,  regardless  of  the  operating 
conditions.  The  differential  equations  describing  such  a system 
contain  terms  in  which  the  constant  time  delay  occurs.  If  these  equa- 
tions are  to  be  studied  by  means  of  an  electronic  analog  computer,  it 
is  necessary  to  have  available  some  device  that  will  provide  the 
delay.  This  device  should  prod.uce  a constant  delay  of  controllable 
magnitude,  and  must  give  no  effect  other  than  that  of  introducing  the 
delay. 

There  are  several  ways  by  which  a constant  time  delay  can  be 
obtained.  For  many  purposes,  however,  it  ia  most  convenient  to 
synthesize  from  the  elements  of  the  analog  computer  itself  a network 
which  will  produce  the  delay.  This  network  is  required  to  transmit 
the  electrical  signal  representing  the  solution  for  the  equation  being 
studied.  Within  the  frequency  band  containing  the  signal,  the  net- 
work should  introduce  constant  known  values  of  bo In  time  delay  and 
attenuation.  Such  a network  is  described  as  an  all-pass  network  with 
its  phase  shift  varying  Linearly  with  the  frequency  of  an  applied 
sinusoidal  signal.  The  linear  relation  between  phase  shift  and 
frequency  should  be  maintained  over  the  necessary  frequency  band. 

Time-delay  networks  of  varying  complexity  can  be  designed.  In 
most  applications  it  is  desirable  to  use  as  few  elements  as  possible 
in  the  network,  and  to  obtain  as  large  a value  as  possible  for  the 
product  of  the  time  delay  x and  bandwidth  to^.  One  design  for  a 
time-delay  network^  has  appeared  recently  in  the  literature.  In  the 


1.  C.  D.  ilorrill,  Trans.  IKS,  EC-3,  45,  (June  195^) 
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following  discussion,  a somewhat  different  approach  to  the  problem  is 
used,  and  severa1  networks  are  analyzed.  The  number  of  elements  and 
the  circuitry  needed  to  synthesize  them  are  shown,  and  the  value  of 
product  <omT  is  obtained.  The  results  of  the  analysis  are  useful 
in  selecting  a delay  network  appropriate  for  a given  application. 

II,  Analysis  of  Networks 
II. 1 All-' pass  networks 

The  networks  considered  here  are  all  of  one  general  type.  They 
are  all-pass  networks,  in  which  the  attenuation  of  a sinusoidal,  signal 
is  independent  of  frequency.  Furthermore,  the  networks  are  adjusted 
so  that  this  attenuation  is  zero,  or  that  the  magnitude  of  the  signal 
is  unchanged  by  the  network.  The  curve  relating  phase  shift  to  fre- 
quency is  essentially  linear  frem  zero  frequency  to  some  maximum 
angular  frequency  For  frequencies  higher  than  this  maximum,  the 

phase  curve  departs  rapidly  from  linearity.  The  time  delay  x is  the 
negative  of  the  slope  of  the  phase  curve 

x **  -dfi/dw  (1) 

where  6 is  the  total  phase  shift.  The  delay  is  constant  so  long 
a; i the  phase  curve  is  linear. 

The  transfer  function  H for  such  a network  is 

H = E2(p)/E1(p)  (2) 

where  Eg  and  are  the  complex  signals  at  the  output  and  input 
terminals,  respectively,  and  p is  the  ooraplsx  frequency  variable. 

In  order  to  fit  the  requirements  listed  previously,  the  magnitude  of 
H should  be  unity  and  its  angle  should  vary  linearly  with  real  fre- 
quency. The  variation  of  H with  p can  be  studied  most  easily 
through  a consideration  of  the  roots  of  the  polynomials  in  p that 
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form  the  numerator  and  denominator  of  the  fraction  i'epresenting  H, 

2 

It  has  uoen  shown  that  these  roots  must  fit  a specified  pattern. 

They  must  occur  either  as  real  quantities  or  as  complex  conjugate 
pairs.  Roots  of  the  denominator,  or  polos,  must  occur  in  the  left 
half  of  the  complex-frequency  plane,  and  roots  of  the  numerator,  or 
zeros,  must  occur  in  the  right  half- plane.  Each  pole  must  have  a 
matching  zero  with  its  real  part  of  opposite  algebraic  sign.  A typical 
situation  is  shown  in  Fig.  1,  where  six  roots  appear.  The  design 
problem  for  the  delay  network  is  essentially  that  of  determining  the 
optimum  location  for  these  roots. 

II. 2 Example  with  six  roots 

The  analysis  of  the  network  with  the  six  roots  of  Fig.  1 is  used 
as  an  example.  The  transfer  function  for  this  network  is 

H = .Cb_- - ia-  ~-4b-8-£  ~ -c2  (3) 

[p  * (a  - jb)J  [p  + (a  + jb)J  fp  + cj 

If  a real  driving  frequency  is  used,  p - ja>,  and  the  angle  of  H is 

0 » tan  '1'(b  - oo)/a  + tan  “(-b  - co)/a  + tan  '1'(-co/c) 

- tan  '1'(-b  + w)/a  - tan  ^(b  + «)/a  - tan  \u)/c)  . (4) 

It  is  convenient  to  normalize  some  of  the  quantities  by  using  the 
definitions 

x a b/a,  y s c/a,  8 s co/a  (5) 

so  thal  Eq.  (4)  becomes 

9 " 2 [tan  ~(x  - 2)  - tan  ^(x  + Q)  - tan~^(ty/y )J „ (6) 

The  procedure  now  is  to  take  successive  derivatives  of  9 with 
respect  to  the  frequency  variable,  Q.  The  first  derivative  is 


2.  H.  W.  Bode,  Network  Analysis  and  Feedback  Amplifier  Design, 
(D.  Van  Nostrand,  New  York,  1945),  p.  239 
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F*g.  1 Location  of  poles  and  zeros  for  6-root  network 
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proportional  to  the  time  do  lay  t.  At  zero  frequency  the  delay  is 

defined  as  t , the  initial  time  delay.  Derivatives  of  even  order 
o 

all  vanish  at  zero  frequency  because  the  phase  curve  is  an  odd  function 
of  frequency.  The  parameters  x and  y are  chosen  so  as  to  make 
the  derivatives  of  third  and  fifth  order  vanish  at  zero  frequency. 

This  insures  that  the  phase  curve  will  be  as  nearly  linear  as  possible 
near  zero  frequency,  and  should  provide  a valid  criterion  for  the 
design  of  the  network. 


The  successive  derivatives  of  Eq.  (6)  are  as  follows. 


p 

It  is  evident  from  Eqs.  (8)  and  (10)  that  at  Q = 0 both  D O = 0 
and  D^9  = 0,  where  0°$  is  defined  as  dn6/’dfln  evaluated  at  0 = 0. 


The  odd-order  derivatives  evaluated  at-  O ■ 0 become 


-2  — 2— r + - 

1 + x2 


1P9  =■  4 


&-.2£U.L 


Ds9  . . Jky:  sA  , f 

[l  * x2]5  y5_ 


The  requirement  that  D^e  * 0 and  D39  = 0 gives 
y3  - -d  ♦ x2)3/2(l  ~ 3x2) 


y5  - -(1  ♦ x2)5/2(1  - lCbc2+  5x4) « 


If  y is  eliminated  from  these  two  equations,  the  result  is 


(1  - lQx2  + 5x4)3  - 4(1  - 3x2)5  . 


The  value  of  x satisfying  this  equation  is  x = 0.953,  end  when 


this  value  is  used  in  either  of  Eqs.  (15)  or  (lo) , the  corresponding 


•value  for  y is  y “ 1.26.  These  numerical  values  serve  to  give  the 


optimum  locations  for  the  roots  of  the  network. 


The  initial  time  delay  is,  from  Eq.  (l). 


t_  --O/a)D0. 


With  the  numerical  values  of  parameters  x arid  y found  above,  the 


time  delay  is  found  from  Eq.  (12)  as 


v “ 3.68/a. 

o 


in  order  to  determine  the  frequency  band  over  which  the  delay 


is  essentially  constant,  it  is  necessary  to  plot  a curve  of  9 as  a 


function  of  £.  Equation  (6),  with  the  numerical  values  of  x and  y 


inserted,  is 


9 - 2[tan”1(,953  - 0)  - tan“1(.953  + Q)  - tan"1(fl/l.26)] . 


This  equation  i?  plotted  in  Fig.  2,  Also  plotted  there  is  the  straight 


line  having  the  slope  found  from  Eq.  (12),  The  curve  cf  Eq.  (20) 
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departs  from  the  straight  line  as  Q becomes  larger.  Purely  arbitrar- 
ily, a difference  of  five  percent  might  be  chosen  as  an  allowable 
departure.  This  departure  occurs  for  Q « 1,3,  Thus,  the  upper 
frequency  limit,  defined  in  this  way,  is  =*  1,3  a,  and  the  product 

a)  t is  a sort  of  figure  of  merit,  w t • 4,8,  An  alternate  form 
mo  'mo 

for  the  figure  of  merit  is  the  value  of  che  phase  shift  when  the 

departure  from  linearity  is  the  arbitrary  five  percent.  Hare,  this  is 

= 277  degrees,  which  is  merely  expressed  in  degrees.  It  is 

the  value  9 would  have  at  o>  if  the  phase  curve  were  truly  linear, 

n 

The  actual  value  of  9 is  five  percent  less  than  this  9^. 

111,3  Other  networks 

All-pass  time-delay  networks,  based  on  this  typo  of  design, 
obviously  can  be  created  with  any  integral  number  of  pairs  of  roots. 


The  example  just  studied  contains  three  pairs. 


nuTubcr 


pairs  is  odd,  one  pair  must  be  located  on  the  real  axis  of  the  p-plane. 
Where  the  number  of  pairs  is  even,  all  of  them  will  be  located  off  the 
axis.  The  complexity  of  the  physical  system  needed  to  produce  the 
required  transfer  function  will  increase  fairly  rapidly  as  the  number 
of  roots  is  increased. 

In  Tables  I-IV  are  given  the  results  of  the  analysis  of  networks 
of  this  sort  having  one,  tv/o,  three,  and  four  pairs  of  roots.  By  the 
time  four  pairs  of  roots  are  used,  the  numerical  work  needed  to  find 
the  optimum  parameters  becomes  quite  tedious.  The  performance  of  this 
8-root  network  is  adequate  to  satisfy  many  applications,  and  more 
complicated  networks  are  not  considered  here. 


f 
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Phase  curves  for  each  of  the  networks  are  plotted  in  Fig,  3 
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Table  Is  2- root  Network 


Pairs  of  roots:  one 

Location: 


(-a  + JO) 


(a  + .10) 


Design  criteria:  D 9 = 0 

Circuit  parameters:  only  a is  specified 


Initial  delny:  = 2/a 


Maximum  frequency:  co^  = 0.4a 


Figure  of  merit:  w - 0,3,  9 = 46  degrees 

° mo  7 m w 


Circuit: 


-<!> 


iO 

‘VI/’  n 


> — 


i 

_.j  i 


x 


Transfer  function:  H = - V 

\P  *r  W 
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Table  II,  1+  - root  Network 


Pairs  cf  roots: 

two 

» 

l 

j 

Location: 

(-afjb)  x 

° (a+jb) 

' 

1 

(-a-jb)  x 

© (a-jb) 

•o  ^ » 

Design  criteria:  D^S  = Ek9  a D*11©  * 0 

for  D36  - 0:  (1  - 3x2)  = 0 


Circuit  parameters:  x * b/’a  » 0,577 

Initial  delay:  t ■ 3/a 

o 

Maximum  frequency:  u = 0,9a 

Figure  of  merit:  wmTo  » 2,7,  = 155  degrees 

Circuit: 


.Mi-  -!»9^ urs****"^ 


• if 

i t- 


Table  III.  6- root  Network 


Pairs  of  roots!  three 
Location! 

(-a+jb)  x 
(-c+jO)  t 


O (a+jb) 


(c+jO) 


(-a-jby  >* 


O(a-jb) 


Design  criteria:  D2G  “ iPo  » D‘+t>  =>  D5G  - DG9  = 0 


for  D3©  = 0:  » 0 


mk 


ii. 


(1  ♦ X2)3  r* 

for  D5Q  - 0:  SiLr_A^-i_S2Ll  + i - o 

(1  + x2)5  y5 

Circuit,  parameters:  x = b/a  = 0.953 

y - c/a  » 1.26 
Initial  delay : tq  = 3.60/a 

Maximum  frequency:  w = 1,3  a 

m 

Figure  of  merit:  w -t  = *+.3.  ? * 277  c agrees 

in  o 'a 

Circuit:  The  circuit  of  Table  I,  designed  with  parameter  c,  followed 

by  the  circuit  of  Table  II,  designed  with  parameters  a and  b. 

Transfer  function:  H - {£=  2aP  * - cj 

1.P  + 2ap  + (a^  + d )]  [p  + c] 


I 


12. 


' 

1 

Table  IV, 

Pairs  of  roots:  four 

Location: 

S-root  Network 

1 

f 

• 

(-c+jd)  * 

° (c+jd) 

i 

(-a+jb)  x 

O (a+jb) 

1 

{ IV  x.  \f  ! 

V-c i-jwy  - | 

A ( - ju  \ 

' 

(_c-jd)  X 

o(c-jd) 

Design 


criteria:  D2S  *■  D39  c =■  D76  = D^#  = 0, 


for  D3©  - 0:  -A.Zj.j-.  + -f-Z.^L-  « 0 

(i  + x2)3  Pa  * i-y 


lor  D-'e  ■“  0: 


1 - 10x2  + 5x*  1 - 10z~  + 5z*  „ 

- ~ — **  o 


(1  + x2)5 
2 

- 21  : 


5m  a 2v5 

y U + z ) 


for  D'6  - 0:  ‘ ^ ; ~ ^ , 0 


r/-»  . _2\7 

jr  \-4-  ~ ; 


(1  + X2)7 

Circuit  parameters:  x = b/a  = 0.34 

y **  c/a  «•-  0*71 
z ■=■  d/c  = 1.5 

Initial  delay:  tq  = 5# 26/a 

Maximum  frequency:  wm  = 1.55a 

Figure  of  merit:  G>mT  = 6,2,  © = 470  degrees 

Circuit:  The  circuit  of  Table  II,  designed  with  parameters  a and  b, 

followed  by  a second  circuit  of  Table  II,  designed  v/ith 
parameters  c and  d. 


[p2  ~ 2ap  + (a2  + b2)]f.p2  - 2cp  + ' -2  + d~)j 
Cp2  + 2ap  + (a2  + b2) j[p2  + 2cp  + (c2  + d2)] 


Transfer  function:  H 


lit 


I 

! 

I 
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III . Synthesis  of  Networks 

The  elements  of  an  electronic  analog  computer  must  now  be  used  to 
synthesize  a network  having  its  roots  determined  by  the  foregoing 
analysis.  Since  tha  computing  amplifiers  are  unilateral  circuit 
elements,  the  synthesis  is  not  difficult  to  carry  out.  One  circuit  i3 
sufficient  to  give  a single  pair  of  roots  on  the  real  axis,  A.  second 
circuit  will  give  two  pairs  of  roots  symmetrically  located  off  the 
axes.  These  two  circuits  can  then  be  combined  to  give  any  arrangement 
of  roots  that  may  bo  desired. 


111,1  One  pair  of  roots 

The  transfer  function  fur  the  time-de'iav  network  with  a single 


pair  of  roots  on  uuc  ruoi  axis  is 


H 


(p  - a 
(p  + a 


(21) 


which  can  be  written  as 

H * 1 - 2a/ (p  * a)  . (22) 

The  tern,  -l/(p  + a),  can  be  obtained  by  applying  resistive  feed- 
back around  an  integrator.  If  this  term,  is  adjusted  in  magnitude  by 
the  factor  ka  and  acidea  to  unity,  the  result  is  the  desired  transfer 
function.  The  corresponding  circuit  is  shown  in  Table  I.  Tha  circuit 
yields  a transfer  function  which  is  the  same  as  Eq,  (21),  but  with  a 
negative  algebraic  sign. 


Ill .2  Two  pairs  of  roots 

The  transfer  function  for  the  network  with  two  pairs  of  roots  is 

2 ~ _ . /_2  . , 2\ 


u „ (p  - a - .ibHp  - a + ,jb)  _ p - 2ap.  + (a  + b ) 
(p  + a > jb)(p  + a -jb)  2 + 2ap  + {2  + b2} 


(23) 
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which  can  be  written  as 

H «*  1 - 4ap/|j)2  + 2ap  + (a2  + b2)J  . (2/..) 

The  term,  + 2ap  + (a2  b2)J , can  be  obtained  by  applying 

both  resistive  feedback  and  feedback  with  one  integration  around  on 
integrator.  If  this  term  is  adjusted  by  the  factor  4a  and  added  to 
’uiity,  the  result  is  the  desired  transfer  function.  The  corresponding 
circuit  is  shown  in  Table  II,  Thi3  circuit  yields  a transfer  function 
which  is  the  same  as  Eq.  (23),  but  with  a negative  algebraic  sign. 

III. 3 Several  pairs  of  roots 

When  several  pairs  of  root3  are  necessary,  it  seems  simplest  to 
ccanbine  the  circuits  of  Tables  I and  II  in  appropriate  fashion.  Thus, 
the  circuits  of  Tables  I arid  II  are  combined  in  Table  III;  two  circuits 
of  Table  II  are  combined  in  Tabls  IV, 

It.  is  possible,  of  course,  to  synthesize  other  arrangements  of 
circuit  elements  that  will  give  the  transfer  functions  vrith  several, 
pairs  of  roots.  The  alternate  configurations  that  have  been  explored 
require  at  least  as  many  circuit  elements  as  the  combinations  suggested 
in  Tables  III  and  IV,  and  appear  to  be  more  difficult  to  adjust. 

Where  a complicated  circuit  is  obtained  as  the  sum  of  several  simpler 
circuits,  each  circuit  can  be  tested  alone  to  make  sure  it  is  operating 
correctly.  Alternate  forms  of  the  complicated  circuit  may  not  be  so 
easily  broken  dovm  for  test  purposes, 

IVS  Practical  Considerations 


The  primary  consideration  in  choosing  the  appropriate  delay  net- 
work for  a particular  application  is  the  total  phase  shift  that  must 
be  obtained  with  a curve  essentially  linear'.  This  is  the  form  of  the 
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figure  of  merit  given  as  6 in  the  tables , A network  vdth  a larger 

1H 

number  of  roots  gives  a larger  value  for  6 . Thus,  if  a large  value 
of  6 is  needed  for  a certain  probler,  a relatively  complicated 
delay  network  is  required. 

Two  identical  networks  might  be  used  in  tandem  to  give  double 
the  delay  of  either  alone,  In  Fig,  4 is  shown  the  phase  curve  for  a 
pair  of  identical  4~root  networks  used  in  this  way.  Also  shown  is 
the  curve  for  a single  3-root  network.  The  initial  delay  for  the  pair 
of  4- root  networks  is  twice  that  for  one  alone,  and  is  larger  than 
that  for  the  8-root  network.  However,  the  departure  from  linearity  is 
greater  than  for  the  single  8- root  network.  The  figure  of  merit  for 
the  pair  of  4-root  networks  is  twice  2,7,  or  5.4,  while  that,  for  the 
single  3-root  network  is  8,2.  This  example  indicates  that  it  is 
advantageous  to  use  a single  network  of  sufficient  complexity  to  give 
the  necessary  total  phase  3hift. 

In  setting  up  tho  components  to  form  a delay  circuit,  it  is  some- 
times helpful  to  test  isolated  portions  of  the  complete  system.  The 
circuit  giving  four  roots  can  be  profitably  broken  down  in  this  way. 

If  the  transmission  path  straight  f ran  tho  input  terminal  of  the  net- 
work to  one  input  terminal  of  the  summing  amplifier  is  removed,  the 
transfer  ^unction  of  the  remaining  elements  is 

G = j4am/jja2  -t  b"  ® to2)  ♦ J2amj,  (25) 

It  is  not  difficult  to  show  that  its  magnitude  is  j G j — 2 for  the 
2 2 2 

driving  frequency  u>  « (a  + b )•  Similarly  its  magnitude  is 
jcU  2^^  for  the  driving  frequencies  u2  - (3a^  + b^)  ± 2a(2a2  + b^)'^. 
This  portion  of  the  network  can  easily  be  checked  to  make  sure  that  J G | 
has  the  proper  value  at  tho  frequencies  given. 


Q 


Fig.  4 Phase  shift  of  two  4-root  networics  and  of  one 
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The  magnitude  of  the  transfer  function  for  any  of  the  complete 
delay  networks  is,  of  course,  unity  at  any  driving  frequency,  and  this 
fact  should  be  checked. 

V.  Example  of  Application 

A nonlinear  differential-difference  equation  requiring  a time 
delay  is 

dX(fc)/dt  =*  A X(t)  - B X(t)  X(t  - t)  (26) 

where  X is  the  dependent  variable,  t is  time,  the  independent 

variable,  A and  B arc  positive  real  constants,  and  a is  a corr- 

3 

stant  time  delay.  This  equation  has  been  studied  on  an  analog 
computer,  using  a 3tep-by-step  process  of  computation.  For  certain 
values  of  product  At,  the  solution  for  the  equation  is  a periodic 
oscillation  with  a waveform  that  is  far  from  sinusoidal.  The  period 
of  the  oscillation  is  at  least  four,  or  more,  times  the  delay  time  t. 
Thus,  a delay  of  sanewhat  less  than  90  degrees  is  needed  at  the  funda- 
mental frequency.  However,  many  harmonic  components  are  present  in 
the  nonsinusoidal  oscillation.  The  phase  shift  for  a harmonic  with 
constant  time  delay  is  the  prcduct  of  the  order  of  the  harmonic  and 
the  phase  shift  of  the  fundamental.  Since  it  appears  that  harmonics  of 
at  least  i,hs  fifth  order  arc  important,  a total  phase  shift  of  at 
least  400  degrees  is  required.  The  use  of  an  8-root  delay  circuit  is 
indicated, 

A convenient  value  for  the  time  delay  with  the  computer  used  here 
is  tq  = 2 seconds.  Then,  using  the  design  equations  of  Table  IV,  the 
necessary  parameters  are 

3,  W.  J.  Cunningham,  A Nonlinear  Differential-Difference  Equation  of 
Growth,  Contract  ilonr-433(00),  ’/ale  University,  May  19V*..  Upper- 
case symbols  are  used  in  •writing  En.  (PA)  +.9  M*lwh 

symbols  used  previously. 


V. 


a -j  5,28/tq  = 2,64  sec 

b = xa  = ,34  k 2,64  * 0,895  see 

c.  = ya  « ,71  * 2,64  “ 1.87  sec  1 

d « zc  ■ 1,5  x 1.87  * 2,81  see  x 

2 2 -2 
+ b ■ 7.72  sec 

+ d’  “ 11 ,41  sec  ^ 

The  configuration  of  elements  needed  to  synthesize  the  delay  circuit 
is  shown  in  Fig,  5«  The  complete  setup  of  the  cccrputor  for  studying 
Eq,  (26)  is  shown  in  Fig,  6,  and  incorporates  the  delay  circuit  of 


~ e 
4 — c» 3 ✓ » 


A family  of  curves  obtained  with  this  system  is  shown  in  Fig.  7. 
The  only  initial  condition  for  these  curves  is  that  X <=-  p/lG  ,t 
t 3 0,  where  £ *=  A/B,  These  curveB  agree  well  with  those  of  Fig.  19 
of  footnote  reference  3»  Those  latter  curves  were  plotted,  with  points 
obtained  from  a step-by-step  integration.  The  curves  of  Fig.  ? are 
plotted  directly  and  continuously  by  the  computer.  The  minimum  value 
of  product  At  needed  to  give  a periodic  solution  is  approximately 


a-.  _ i rtc  muj  - -•  » M 
Ml  • A J,  1X0  xo  uumo  n 


what  larger  than  that  estimated  in  the  reference , 


The  solution  just  overshoots  its  final  value,  so  oscillation  begins, 
if  At  = 0,48,  which  is  close  to  the  previous  estimate. 

In  Fig.  8 is  shown  the  solution  for  X(t)  with  At  “ 1.3,  and 
also  this  same  function  as  retarded  by  the  delay  network,  Except  for 
the  initial  transient,  the  delayed  curve  is  quite  accurately  the  original 
curve,  merely  shifted  the  proper  amount  on  the  time  scale.  If  ohe 
product  At  is  much  larger  than  about  Ac  * 2,  the  oscillation 
becomes  even  more  noneinusoidal.  The  higher  harmonics  then  exceed 


fhp  mA-sHrrunn  frequency  for  the  delay  network,  and  the  delayed  signal 
begins  to  differ  in  shape  from  the  signal  with  no  delay. 
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Fig,  6 Cofnputer  aelup  for  Eq«  (26) 
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The  initial  transient  on  the  delayed  curve  of  Fig.  8 ecmes  about 
because  no  initial  conditions  were  applied  to  the  integrators  of  the 
delay  circuit.  If  enough  information  is  known  about  the  solution  oi 
the  equation  for  the  interval  -x  <£  t < 0,  presumably  initial  condi- 
tions needed  in  the  delay  circuit  could  be  predicted  and  inserted  in 
such  a way  as  to  avoid  this  initial  transient.  In  practice,  however, 
problems  involving  differential-difference  equations  of  this  sort  appear 
to  have  the  inherent  difficulty  of  not  furnishing  sufficient  informa- 
tion to  start  the  solution  properly.  As  has  been  observed  previously, 
this  kind  of  equation  is  equivalent  to  a differential  equation  of 
infi  nitiC  orde  r,  so  that  an  infinity  of  initial  conditions  arc  needed, 
with  a finite  delay  network,  such  as  those  considered  here,  only  a 
finite  number  of  conditions  can  be  used.  Evon  this  finite  number  ’dll 
not  be  known  in  the  general  problem. 

The  damping  effective  in  the  delay  circuit  is  sufficient  so  that 
if  there  is  a transient  set  up  within  it,  decay  is  fairly  complete 
within  the  interval  of  one  delay  time,  t. 

Initial  conditions  needed  to  eliminate  the  initial  trar-  'v  .1  in 
the  -3-root  delay  network  can  bo  found  as  follows.  The  network  with 
its  four  integrators  is  shown  in  Fig,  9.  Initial  conditions  applied  to 
the  integrators  are  identified  as  IC  ; these  are  the  negatives  of 

fit 

the  sigiials  which  must  exist  at  the  output  terminals  of  the  integrators 

at  zero  time.  Constant  multiplying  factors  for  certain  of  the  computing 

amplifiers  are  identified  as  N , Signals  at  various  points  in  the 

n 

systems  are  called  = e^(t). 

Through  a consideration  of  the  transfer  functions  for  the  parts 
of  the  network,  the  following  equations  can  be  written. 


Fig.  9 Insertion  of  linitial  conditions  in  8-root  network 
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P (4/N,)ap  “|  j 

- = -e.(0)  » -5 2~ — ~2~  | j 

1 2 j_p2  + 2ap  + (a2  + b2)J  1 |t-0 

N2e2(0)  j"  -(W/Nj.Ja  1 j 

2 P L»2  + 2ap  + (a^  + b2)J  x (t-0 


(Vf’o)cp 

I0o  * -e,  (0)  “ —x “~5 — -7T- 

^ ^ Lp  ♦ 2cp  *»•  (c  + d ) 


I o 9 n 

p ' — 2ao  + (a  * ♦ d ) 


Mt)! 


5 0 | —i  v "i  i 

+ 2&p  + (a^  + d ) (t=0 


„ v»(°)  r 

**  p |_p2  + 2cp  + (c2  + d") 

L p2  - gap  ♦ (a3  » b^)  j 

1.  p2*2ap<  (,S  >b2)j  ’'“'It-Q 


Enough  information  is  known  about  Eq.  (26;  to  detoruxne  th^se 
initial  conditions.  Near  X «»  0,  approximately  X = X exp (At), 
where  X =>  X_  at  t =*  0,  For  the  curves  cf  Fig.  8,  A - 0*9  and 
Xq  « <3/10.  T.tu3,  for  the  3etup  of  Figs.  5 and  6,  and  the  curves  cf 
Fig,  8,  numerical  data  needed  for  Eqs,  (27- ’30)  are  as  follows: 


a = 2*64 
b » 0,895 
c = 1,87 


d • 2.81 


H1  “ '» 
N2  = !0 


N,  » 10 
4 


^(t)  « (p/l0)exp(0.9t) 


When  these  data  are  inserted  into  the  equations  and  the  necessary 
operations  carried  out,  the  initial  conditions  for  the.  integrators 


are  found  as 


1^  = 0.177  (3/10) 
IC?  - -1.97  (8/10) 


IC3  « -0,0123  (8/10) 
IC.  » 0.137  (p/10) 


When  these  conditions  are  applied,  the  delayed  solution  starts  off 
smoothly,  without  an  initial  transient  such  as  that  of  Fig.  8. 
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